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We present a dimensionally deconstructed model of an s-wave holographic superconductor. The 
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I. INTRODUCTION 



High-T c superconductors continue to fascinate physicists, both for their technological 
prospects and for their unconventional physical properties. Discovered in 1986 jlj], high-T c 
superconductors have superconducting transition temperatures T c higher than the upper 
limit of around 30 K suggested by the BCS theory {2], and the mechanism of superconduc- 
tivity in these materials remains uncertain. High-T c superconductors are Type II, allowing 
magnetic fields to penetrate in conjunction with Abrikosov vortex currents which maintain 
localized magnetic fluxes. A number of features of high-temperature superconductors re- 
main poorly understood, including the existence of a pseudogap phase in which a gap in 
the excitation spectrum opens up at temperatures above the superconducting transition 
and a thermoelectric Nernst effect that occurs in both the superconducting and pseudogap 
phases {4]. A few additional classes of exotic high-T c superconductors have recently been 
discovered, including the iron pnictides [5[, and research into their properties is ongoing. 

In the absence of a clear theoretical understanding of high-T c superconductors and other 
systems with strongly correlated electrons, it is useful to explore simplified models which 
describe similar phenomenology. It has long been understood that certain relativistic field 
theories describe a variety of phenomena typical of superconductors. In the Abelian Higgs 
model, a charged scalar field condenses while the photon becomes massive, in analogy to 
the condensation of Cooper pairs and the consequent expulsion of magnetic fields from su- 
perconducting materials. The Abelian Higgs model also has Nielsen- Olesen vortex solutions 
which maintain localized magnetic flux tubes 6[, much like the Abrikosov vortices in Type 
II superconductors. 

One approach that has stimulated much effort especially from the string theory commu- 
nity is the use of the AdS/CFT correspondence which relates certain strongly coupled 
field theories to gravitational systems in higher- dimensional spacetime backgrounds. Due 
to the extra-dimensional nature of these models they are referred to as holographic. There 
is a mapping between observables in the superconducting system and fields in the corre- 
sponding extra- dimensional model. Progress has been made in understanding properties of 
holographic models of superconductors and other condensed matter systems, including s, p , 
and c?-wave superconductivity sllQ], the Nernst effect in the quantum critical regime 
strange metallic behavior 12|, enhancement of the superconducting gap compared to BCS 



2 



theory 13] . and Homes's empirical law for T c as a function of low-temperature superfluid 
density and normal-phase resistivity [l^, flol. In some models of holographic superconduc- 



tors there is no hard zero-temperature gap 
condensed matter systems, see Refs. [l3, LL2 



16 L For some reviews on holographic models of 
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In parallel to the development of the AdS/CFT correspondence in string theory, extra- 
dimensional model building has led to the discovery of new paradigms for gravitational 



physics 
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20], electroweak symmetry breaking 21], grand unification 22], and other pos- 



sible new physics. Notably, in order to define these otherwise nonrenormalizable extra- 
dimensional gauge theories, a procedure was developed for dimensionally "deconstructing" 
extra-dimensional theories into renormalizable field theories defined in the natural (lower) 



dimension of the system of interest 



23] . At long distances, a deconstructed model is ef- 



fectively described by the corresponding extra-dimensional theory, latticized in the extra 
dimension(s). In some models the number of effective lattice sites can be taken as small 
as just two or three while retaining the desirable features of the extra-dimensional theory. 
Deconstruction provides a systematic procedure for finding relatively simple, weakly cou- 
pled models of interesting physical systems, in the number of dimensions that naturally 



characterizes the system of interest 



24). 



In this paper we present and analyze a deconstructed model of the simplest holographic 
superconductor, the Abelian Higgs model in 3+1 dimensional 1 Anti-de Sitter-Schwarzschild 
spacetime (AdS4-Schwarzschild), which is intended to model s-wave superconductors in 2+1 
dimensions |9|. Superconductivity in the cuprates is along two-dimensional CuC>2 planes, 
so the hope is that a model with two spatial dimensions will capture some of the relevant 
physics in these systems. The deconstructed model is a 2+1 dimensional model in which 
Maxwell's electrodynamics is replicated a number of times. Fields charged under pairs of the 
electromagnetic U(l) gauge groups condense, giving rise to the latticized extra-dimensional 
structure of the model. The temperature-dependent couplings of the Maxwell fields are 
rotationally invariant but not Lorentz invariant. Additional charged fields representing 
Cooper pair operators condense in the superconducting phase, and the combined effects of 
the multiple condensates determine the properties of the superconductor. The electron pair 



1 The +1 in the spacetime dimensionality always refers to time, not the extra spatial dimension of the 
holographic model. 
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wavefunctions in high-T c superconductors are thought to have more complicated symmetry 
than the s-wave described by scalar Cooper pair fields, so the present model is not expected 
to capture features sensitive to the symmetry of the wavefunctions. 

The deconstructed model helps to elucidate certain aspects of holographic superconduc- 
tivity and suggests the relevant effective degrees of freedom in physical realizations of these 
superconductors. One generic feature of this class of models is the existence of hidden lo- 



cal symmetries 



25| and corresponding neutral excitations, much like the vector mesons of 



Quantum Chromodynamics. Such excitations in the superconductor may be interpreted as 
excitons 2 , electron-hole bound states which can be created when a photon is absorbe d by 
certain materials 



26l | . In analogy to a deconstructed holographic model of hadrons 



we derive discretized AdS / CFT relations between bulk and boundary observables. We ana- 
lyze properties of the superconducting transition and the frequency- dependent conductivity 
as the number of extra-dimensional lattice sites, and correspondingly the number of pair 
condensates and exciton fields, is reduced. 

In Sec. HIl we review the Abelian Higgs model in the AdS4-Schwarzschild spacetime 
as a model of superconductivity in two spatial dimensions. In Sec. IIIIl we describe the 
deconstructed model and explain the calculation of observables of interest. In Sec. IIVt we 
present numerical results. We conclude with a discussion and suggestions for future research 
in Sec. El 



II. CONTINUUM MODEL 

In this section we review the holographic model that will be deconstructed in Sec. IIIIl 
Much of the content of this section is a summary of results contained elsewhere, for ex- 



ample in Ref. js], though we generalize certain results for the sake of comparison with the 
deconstructed model. 

The starting point is a vacuum solution to Einstein's equations with negative cosmological 
constant in 3+1 dimensions, namely the AdS^Schwarzschild solution. The spacetime is 



2 The neutral, spin-1 excitations might more appropriately be identified with polaritons, if there were a 
dynamical photon in the deconstructed theory that mixed strongly with these excitons. 
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described by a metric of the form 



ds 2 = F(r)dt 2 - r 2 (dx 2 + dy 2 ) - -^-c/r 2 , (2.1) 

F(r) 



where 



7 ' 



F (r) = T2 {l-^). (2.2) 



The constant L is the Anti-de Sitter scale, which we will often set to 1. We refer to coordi- 
nates in which the metric takes the form of Eq. (12. ip as r coordinates, and the coordinate r 
runs from r# to oo. The 2+1 dimensional surface at r = oo is referred to as the ultraviolet 
(UV) boundary, and the surface at r = r# is the horizon. The AdS/CFT correspondence 
suggests the identification of the temperature of the material with the Hawking temperature 
of the AdS black hole, namely 

T = 3r H /(47rL 2 ). (2.3) 

With nonvanishing charge density the relevant background spacetime is instead the AdS- 
Reissner-Nordstrom spacetime, a charged solution to the Einstein-Maxwell equations with 
a negative cosmo logical constant. However, in this paper we neglect the backreaction of the 
charge density on the geometry. In this approximation the spacetime metric continues to 
take the form of Eqs. ( 12 .ip and ( 12. 2f) . even when the solutions of interest have nonvanishing 
charge density. 

The holographic model contains a U(l) gauge field corresponding to electromagnetism, 
and a charged scalar field corresponding to the Cooper pairs. The model is translationally 
invariant in 2 spatial dimensions, and is therefore not expected to reproduce phenomena 
that are sensitive to the atomic lattice. 

The action for the model is given by 

S = J d 4 x^{~F MN F MN + \(d M -iA M )i>\ 2 -m 2 \i>\ 2 y (2.4) 

where index contractions are with the 4D metric gMN, and g = |det Qmn\- 
The action in r coordinates has the form 

S=f Sxdrr 2 l^rAFoa) 2 - V 4 (F afc ) 2 + i(F 0r ) 2 - ^ 2 {F ar ) 2 

1 , (2-5) 

+ -— |(<9o - tAM 2 - r- 2 \(d a - iAM 2 - F(r)\(d r - iA r )^\ 2 - m 2 \^\ 2 . 
F{r) J 

Here Fmn = 9mAn — d^AM is the field strength of the U(l) gauge field, and the charge 

of the field ip is normalized to 1. Our convention is that Latin indices a, b, etc., represent 



spatial components x or y. We have explicitly shown the metric factors in Eq. (12. 5p . so 
remaining index contractions are taken with the Kronecker delta 5mn- According to the 
AdS/CFT correspondence the mass m is related to the scaling dimension of the operator 
related to ip, and for definiteness we will assume m 2 = —2/L 2 as in Ref. [9|, corresponding 
to a Cooper pair operator of dimension 2 (or 1, as discussed below). 

Defining = Aq as in Ref. 9J, and considering solutions with ip = if)(r), — 0(r) and 
A a = 0, the coupled equations of motion for ip and in the gauge where A r = are 

„ fF'ir) 2\ . 2 , m 2 , 

Near the UV boundary, the solutions behave as [9|, 

^ = ^— + ^- + ..., 2.8 

= A* — -H ■ (2-9) 

r 

The time component of the gauge field couples to the charge density, so the AdS/CFT 
correspondence identifies \x with the electric chemical potential and p with the charge density. 
For the Cooper pair operator there are two choices: either acts as a source and ip^ 
corresponds to the expectation value of the Cooper pair operator, or vice versa. This choice 
determines the scaling dimension of the Cooper pair operator according to the AdS/CFT 
correspondence. For definiteness we will choose to take ip^ as the source of the Cooper pair 
operator O2, which then has mass dimension two; the other choice would correspond to an 
operator of dimension 1. The Cooper pair operator should only be turned on dynamically, 
so vanishing of its source becomes a UV boundary condition, tp^ = 0. Additional boundary 
conditions are <f>(rn) = for regularity of the solution and 0(oo) = //. The equations of 
motion enforce the condition F'(r)ip'{r) = m 2 ^(r) at the horizon if ip(r) is finite, so the 
additional boundary condition on ip at the horizon is replaced with a regularity condition. 

In order to study the response of our system to an oscillating electric field, we consider 
solutions for which A x oscillates in time, A y = A r = and Aq behaves as in Eq. (12. 9ft . 
Assuming the form 

A x (t,r) = e~ iut A{r), (2.10) 



6 



the bulk equation of motion for A x has the form, 

• 2 d 
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F(r) dr 
For large r the solutions for A x have the form 



A(r) - — (F(r)A'(r)) + 2A(r)\ip{r)\' 



0. 



(2.H) 



A x = Af+A x ^/r + 



(2.12) 



According to the AdS/CFT correspondence, A x ^ acts as the source for the current j x and 
hence corresponds to a background electric field E x = dtA x °\ The AdS/CFT correspondence 
also suggests an identification of the normalizable component Ax with (minus) the current 



j x in the given background state 
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]. In Sec. IHIt we explain how this identification arises 
in the deconstructed model. Hence, we can write 3 



E x = d t A x \ 



-iu)A x \ 



j r i i x r^oo . 



f = -4 1 ) = r 2 d r A x 

The conductivity a is therefore determined in the holographic model by 

r 2 A'{r) 



(2.13) 



a 



Err 



-iuA(r) 



(2.14) 



The Normal Phase 



At temperatures greater than T c the charged field i/j vanishes, in which case solutions of 
the following first order equations also solve Eq. (12. Ill) : 



The solutions are 



^4±( r ) = ex P |±2W J 



exp 



1UJ 



6r H 



A' ± {t) 



dr 

2V / 3tan" 1 



/Co' 



A + (r). 



(2.15) 



2r + ru 
V3r H 



log 



r - r H 



r 2 + tth + r H 



(2.16) 



and the generic solution to Eq. (12. lip is a linear combination A(r) = c + A + (r) + c_v4„(r) 
with constants c + and c_. The choice of branch of the inverse tangent and the logarithm in 



3 The sign differences in these expressions compared with Ref. [9( are due to differences in convention reiated 
to the signature of the metric. We follow the conventions of Jackson's Classical Electrodynamics [2F 
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Eq. (12.16P determines an w-dependent constant multiplying each of the two solutions, which 
can be absorbed in the coefficients c+ and c_. The solution A_(r) describes a wave flow- 
ing into the black hole horizon, and A + (r) describes an outgoing wave. Son and Starinets 
have advocated the choice of ingoing-wave solution as a boundary condition at the horizon, 
based on the requirement of causality of correlators calculated by the AdS / CFT correspon- 



dence (30|. For now, we explore the consequences of the generic solution, for guidance as 
to what to expect from the deconstructed model. A natural choice for the boundary condi- 
tion in the discretized theory is more ambiguous because causality may be imposed at the 
level of the lower-dimensional theory, without regards to the effective higher-dimensional 
description. 

The complex conductivity as a function of frequency distinguishes the normal and su- 
perconducting phases. A hallmark of perfect conductivity is a delta function in the real 
part of the conductivity at zero frequency which, by the Kramers-Kronig relation for the 
conductivity, is tantamount to a zero- frequency pole in the imaginary part. The Kramers- 
Kronig relation follows from the assumption that correlations are causal. As a consequence, 
the conductivity ct(lj), which by definition is the Fourier transform of the response function 
relating the current to a background electric field, is analytic in the upper half plane. The 
Kramers-Kronig relation is then the statement of Cauchy's integral theorem for the following 
integral: 



dJ = V / v 1 dJ - in a(cj) = 0, (2.17) 

where the contour C spans the real axis above the pole at u/ = uj and closes in the upper 
half plane, and V represents the principal value of the integral. The conductivity <j(u) is 
assumed to fall off in the upper half plane faster than l/|w|, so that that integral over the 
contour at infinity vanishes. As a consequence of Eq. (12.171) . if the real part of c(uj) has a 
delta-function at u = then the imaginary part has a pole at u — 0, and vice versa. 

It follows from Eqs. (I2.14p and (12.151) that when ip vanishes, the zero-frequency 
pole in Im a is generically absent; the holographic description indeed describes a non- 
superconducting phase when the condensate vanishes, for a generic choice of boundary con- 
dition at the horizon. Note that for both the ingoing-wave and outgoing-wave solutions, the 
normal-phase conductivity is independent of u>: a = +1 for the ingoing wave, and a = — 1 
for the outgoing wave. This result provides a phenomenological motivation for the ingoing- 
wave boundary conditions: the outgoing-wave solution would describe an unusual situation 
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in which the current produced by an electric field points in the direction opposite to the 
electric field. Note also that the conductivity in this model does not indicate the existence of 
a pseudogap phase, and does not fall off at large u due to relaxation as in ordinary materials. 

For the generic solution for A x with vanishing if) and with the identification r# = 4irT/3, 
Eq. (I2.14p gives the normal-phase conductivity: 

c_exp(^)-c + exp(^) 

a(co) = ) ( ) f. (2.18) 

c _ exp +c+ exp(^ 

The real part of the conductivity is then, 

Rea(c)= |C -' 2 ;'!£ (2 ' 19) 

|c_p + \c + \ z + c_c?j_e 4T + c + c_e 4T 

which oscillates as a function olu/T around a = 1 for constant |c + | -C |c_|. 



B. The Superconducting Phase 

For small enough T, the coupling between <fi and i/j leads to an instability which turns 
on the Cooper pair condensate (O2). The condensate depends on the chemical potential [i 
and the charge density p determined by the solution to the coupled equations of motion. As 
argued in Ref. 9J], the critical temperature scales as p 1 / 2 . Near T c the condensate behaves 
as (O2) oc (1 — T/T c ) 1//2 . The complex conductivity for generic T < T c has a delta function 
at oj = and a gap u g , as expected for typical superconductors, but with u g /T c pa 8 {{]], 
significantly larger than the BCS prediction uj 9 /T c 3.5. The phenomenology of the decon- 
structed model has similarities to that of the continuum theory, even with only a few lattice 
sites, as we will see in Sec. [TV] 

III. DECONSTRUCTION 

In the deconstructed model we will define coordinates differently, replacing r with a new 
coordinate z = 1/r. In z coordinates, the extra dimension corresponds to a finite coordinate 
interval. The metric in z coordinates has the form 



ds 2 = - 2 

z l 



f(z)dt 2 — -jj-rdz 2 — (dx 2 + dy 2 
J \ z ) 



(3.1) 
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where 

1 



^H^li-prh (3-2) 



where z runs from to Zh- The surface z = is the ultraviolet boundary, and z = zh is the 
horizon. One should keep in mind that results are coordinate independent in the continuum 
model, but away from the continuum limit, deconstructed models depend on the choice of 
coordinates in which their extra-dimensional parent model is defined. 
In z coordinates, the action in the continuum theory takes the form 

S = J d'x^-Fl + ±- f Fl - f -F 2 a - \Fl + -L|^ - 



-L\d z 4, - iA z ^\ 2 - \\d a iP - iA a 4,\ 2 - - A m 2 \^\ 2 \ , 

z A z A 



where repeated a or b indices are summed over the spatial coordinates x and y. As before, 
we use the notation of Ref. js], where we define = Aq. We replace the z axis by a lattice 



of N points, 



e+(j-l)a forj = l...iV-l 

zj = { (3.4) 
e + (N - 2) a + a H for j = N , 



where z^ = zh and e is a UV cutoff. The lattice spacing aj between the j and j + 1 th lattice 

sites is therefore a for j — 1 ... N — 2 and an for j — N — 1. We give ourselves the freedom 

to vary an away from a to test the sensitivity of our results to the density of sites nearest 

to the horizon. The discretized action is given by 

N-i f , 

S = a i J d " x {^ d ° A *3 ~ #j? + ^( F oa)? - f (d a A zj - A! aj f 



f. 

-j{Fab)'j + -hr\doi>j - i<t>jipj\ 2 - 41^- - iA zj ^ 2 

<i Zj Jj Zj 



-\\d a ipj - iA aj ipj\ 2 - \m 2 \ipj\ 2 1 

Z, Z, ) 



(3.5) 



J 3 

where fj = f(zj). We define derivatives in the discretized action such that 



= (0i+i - <t>j)l a 3 > ( 3 - 6 ) 

and similarly for the other fields in the theory. We now construct a 3D theory that reproduces 
Eq. (13. 5p . up to differences at the j = 1 and N boundaries. These differences will be required 
for a proper holographic interpretation of the deconstructed theory. 
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j=1 J=2 j=3 j=N 

UV HORIZON /IR 

FIG. 1: Moose diagram representation of the deconstructed theory. 

Consider a 3D theory with N global U(l) symmetries of which N — 1 are gauged. We 
label the N — 1 gauge factors by the index j = 2 . . . N. We introduce a set of bi-fundamental 
complex scalar fields Sj, j = 1 . . . N — 1, where Ej transforms under the j th and (j + l) th 
U(l) factors, with charges —1 and +1, respectively. In addition, we introduce N complex 
scalar fields ipj, j = 1 . . . N, where ipj transforms under the j th U(l) factor, with charge +1. 
The particle content and charge assignments of this theory are conveniently summarized by 
the moose diagram shown in Fig. [U The solid circles represent U(l) gauge groups, while 
each line connecting to a given circle represents a complex scalar field that transforms under 
that group. The global U(l) symmetry is represented by the dotted circle at the left end of 
the moose. 

We assume that the bi-fundamental fields develop vacuum expectation values (vevs) 
(Ej) = vj. Although the first link field Ei is charged under the global U(l) factor, the 
global symmetry nevertheless remains unbroken if the ipj have vanishing expectation values. 
This can be seen by noting that the effect of a global phase rotation on (Ei) can be undone 
by a (constant and spatially uniform) U(l) gauge transformation at the j = 2 site. The effect 
of this phase rotation on (E 2 ) can then be undone by a U(l) gauge transformation at the 
j = 3 site, and so on. In this way, one sees that the global symmetry of the first site remains 
an invariance of the vacuum, even when the link field vevs assume a nonvanishing profile. 
We therefore identify the current associated with the global symmetry of the first site as the 
QED current j M . As has been discussed in the literature 9|, even though the U(l) symmetry 
associated with QED is not gauged, the model still describes aspects of superconductivity to 
the extent that a treatment of electromagnetic fields as external backgrounds is appropriate. 
In the deconstructed model we could simply add a kinetic term for the associated U(l) gauge 
field to describe a fully dynamical electromagnetism, but in order to match the AdS/CFT 
description we will not do that here. 
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In order to compute expectation values of j^, we will turn on a non-propagating gauge 
field at the first site, A±, that couples to this current. We will also assume that the gauge field 
at the last site A^ N is non-propagating and assumes a specified background configuration; this 
will allow us to impose incoming-wave boundary conditions near the horizon. The boundary 
fields ip\ and ipN are also assumed to be non-propagating below, but only to simplify the 
discussion. The case in which all the ijjj are dynamical is discussed in an appendix. 

The Lagrangian for the 3D theory in Fig. [1] is of the form 

JV-l r ^ ~i N-l 

4 1 



N-l r i N-l 



j=2 



3=1 



where Vj determines the scalar potential at the j th site. The coefficients Zj and the 3D 
metric g^ u vary from site to site. Although we use the term "metric" to refer to g^ u , the 
g^ v simply encode the Lorentz- violating couplings of the theory. The covariant derivative is 
given by 

= <9%- + iAfij - iAt+Jlj , (3.8) 

and 

D^ipj = d^j - iA^j . (3.9) 

Notice the absence of and Aj kinetic terms for j = 1 and j = N in Eq. (13. 7p . We identify 
the non-propagating fields that appear at the ends of the moose with fixed backgrounds in 
the ultraviolet (UV) and the infrared (IR): 

At = A» uv (xn, ik = ifovV), 
A N = Al R (x") , fa = ^m{x v ) . 

We can now compare Eq. (13 .7p to the Lagrangian obtained by latticizing the z coordinate 
of the continuum theory, Eq. (13. 5p . Up to differences that vanish in the a — > limit, one 
finds after some algebra that the latticized theory is recovered if one chooses 

2?° = T' 9f = S ab , Zj = \ and Vj = -M , (3.11) 
the scalar potential 

2\vjilj j+1 - tjjjEjl 2 + m^iptf/z] for j = 2...N -1 
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and the identification 



Vj exp[icijA Z j] ~ Vj(l + icijAzj] 



(3.13) 



The potential is not the most general one consistent with the symmetries of the theory; as 
in any deconstructed model, the form of the action is dictated by the requirement that it 
reproduce the latticized action obtained from the continuum theory. For example, Eq. ( 13 . 1 2 j) 
is fine-tuned so that the (ip'j) 2 term in Eq. (13 . 5 p is reproduced when is set equal to its 
vev. For definiteness, we fix the mass parameter m 2 to its holographically motivated value 
m? = —2. One could allow m to deviate from this choice in more general theories, but we 
will not consider that possibility here. The temperature dependence of terms in the action 
is inferred from the form of fj, which depends on temperature via zh = 3/(47rT). One could 
also imagine more general moose models in which the temperature dependence of the fj 
is determined directly from the microscopic properties of the system. Here we will strictly 
consider the fj that follow from discretizing the continuum holographic theory. 

In what follows, we work in unitary gauge, corresponding to the gauge A z = in the 
four-dimensional theory. As in Eq. (I3.13p . we ignore the physical fluctuations of the link 
fields about their vevs. Then, the Lagrangian of the moose model may be written 

N-l 



C 



N-l 
N-l 



'|2 



3=2 



3=2 



\doipj - i<Pj^j\ 2 2 \ d ^j ~ iA aj ipj\ 2 ^m 2 !^ 



Z 3 



(3.14) 



One may now derive the equations of motion for the dynamical fields (j = 2. . . N — 1), 
assuming the same ansatz applied in the continuum theory. For <f>j that are time-independent 
and spatially constant, the equations of motion for the dynamical fields are given by 

2 



ft: 



o 



(3.15) 



where </>'• 



- 2<f>j + <f>j-i)/a 2 , for j 



j J 3 

f ™ - - 2 . . . N 



2 and 



[{4>n - 4>N-i)/aH 



m-i 



l> N-2)/ a \/ a H- F° r the same ansatz, the equations of motion for the x/jj are given by 



z l fj- 



m 



"3 \ A j-1 Jj J Jj ^jJj 



l/jj = 



(3.16) 



where -0" is defined analogously to cj)" and = (if)j+i — ipj)/dj. When <pj is nonvanishing, 
corresponding to a nonvanishing chemical potential and charge density, the coupling between 
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<pj and ipj in Eq. (13. 16[) acts as a negative squared mass term for each Cooper pair field ipj, 
creating an instability that is enhanced for zj near zh by the temperature-dependent factor 
of l/fl 

Finally, we will require the equation of motion for the spatial components of the gauge 
field, assuming spatially constant fields and the time dependence 

A aj {x^=A aj e-^. (3.17) 

Again for j = 2 . . . N — 1, one finds 



A - + ^V~f)^ + \Jl' 2 ^) A ^°- (318) 

where A" and Aj are defined analogously to i/j'J and ip' y The exciton fields (A a )j are excited 
collectively by the external electromagnetic field (A a )i. 

The horizon boundary condition on Eqs. (13. 15ft and ( 13. 16ft that correspond to those of 
the continuum theory are 

2 

(p N = and il)' N _ x = - — ip N . (3.19) 

6Zn 

If the ipN field were dynamical at the last site, it turns out that the same boundary condition 
on ip would follow from the ipx equation of motion in the continuum limit, as we discuss in 
the Appendix. The boundary condition on the spatial components of An can be chosen to 
reproduce the incoming- wave boundary conditions in the continuum limit, as we discuss in 
more detail below. 

In the continuum holographic theory, physical quantities of interest are related to the 
values and derivatives of the fields at the UV boundary. In the deconstructed theory, we 
find that similar relations apply to the expectation value of the charge density p and the 
current density j a , which are identified via their coupling to the background field Ajj V : 

, N-l 



f = i- 



s 



This reduces at tree-level to 



in |/ n vA Pi ) i e j • (3- 2 °) 



r = - j J—S A [A uv \, (3.21) 

where the dependence on Ajjv arises since the fields are evaluated on their classical equations 
of motion. In the linearized theory, one finds algebraically that S c i evaluates to a sum of 
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surface terms 

S d = J d 3 x [C uv + C IR ] , (3.22) 
where the terms involving the exciton fields are given by 

C uv = M<h - <k) + ^/l A al {A a2 - A al ) (3.23) 

and 

ClR = — /jV-l AiaKAiN — ^aJV-l) , (3.24) 

using the boundary condition (fix = 0. Before evaluating Eq. ( 13.21ft . we first arrange that the 
IR surface term vanishes. This is consistent with the approach in the continuum holographic 
theory, where the IR surface term is also present, but is simply discarded To eliminate 
Eq. (13.24ft . we can add an ad hoc term to the Lagrangian such that the desired IR boundary 
condition is consistent with vanishing IR surface term. Near the horizon, the ingoing- 
wave solution satisfies Eq. ( 12.15ft . which can be imposed as a boundary condition in the 
deconstructed model in a number of ways. We choose to approximate the ingoing-wave 
boundary condition as 

A'aN-n-l = "7 AiTV-n (3.25) 

JN-n-1 

for some n, with n <C N in the continuum limit iV — > oo. If n is taken too small, the factor 
of 1/f in Eq. (13.25ft is large and magnifies the discrepancy between the discretized and 
continuum boundary conditions. Integrating from j — n towards the horizon and solving 
the cu 2 -dependent equations of motion, one obtains a relation between A aN and A aN _i that 
generally leads to a nonvanishing surface term. However, if we add a new term to the 
Lagrangian of the form 

C H = t[c{u 2 )A aN -A aN „ 1 ]\ (3.26) 

then the new IR surface term, evaluated on the solution to the equations of motion, be- 
comes a function of c(u 2 ). This function can be chosen so that the surface term vanishes. 
Eq. ( 13.26ft can be expressed as a function of the gauge-invariant field strength tensors F xz n-i 
and Fo X N-i, and the function c(u 2 ) may be interpreted as a function of — d 2 , which ulti- 
mately acts on the non-dynamical field A a ^. The function c(u 2 ) may be replaced by a 
polynomial approximation, valid at least over some finite range in u. In this way, our dis- 
crete approximation to the ingoing- wave boundary condition, Eq. (13. 25 p . can be imposed n 
sites away from the IR boundary, while the IR surface term is arranged to vanish. 
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The expectation value of the current thus depends only on the UV surface terms, as in 
the continuum holographic theory. Each j = 2 field in Eq. (I3.23P can be related to the 
corresponding boundary field by the discrete version of a bulk-to-boundary propagator. For 
example, we could write <pj = Bjifa, for some bulk-to-boundary propagator Bji. Then the 
charge density is found by computing 



5 



—^(Bn-l) 



1 \ (B 21 - 1) . (3.27) 



a 

The quantity 0i(i?2i — is nothing more than what we would find by numerically solving 
the <p equation of motion and evaluating 4>'\ as we defined it previously. Hence, we conclude 

P = , (3-28) 

in agreement with the holographic prescription for computing the charge density in the 
continuum theory. By the same reasoning, the current density following from Eq. ( 13 .23 j) is 
given by 

f = -fiA' al , (3.29) 

which also agrees with the usual holographic prescription in the continuum, where f± — > 1. 
Since j° is proportional to the charge density operator, we identify <f>\ with the chemical 
potential 

H = (j>i. (3.30) 

The transition to the superconducting phase occurs when the ipj develop a nonvanishing 
profile. In the continuum limit, ijj(z) takes the form ip(z) = ip( z 2 near the UV 

boundary, so that the existence of the condensate is determined by nonvanishing first or 
second derivatives: 

^ = i/;'(z = 0)^0 or ^ (2) = ^"(z = 0) ^ , (3.31) 

In the discretized theory, we use the analogous expression ipj = tp^Zj + ip^ 2 'Zj evaluated 
at the first two lattice sites {z\ = e and z 2 = e + a) to define the operators ip^ 1 ' and 

^e 2 -^(e + a) 2 (2) _ ^ 2 e-^(e + a) 

^ ~ e 2 (e + a)-e(e + a) 2 * " e(e + a) 2 - e 2 (e + a) ' {dJ2} 

These definitions have the appropriate continuum limit. As in Sec. Ull we restrict ourselves to 
the case where ip^ = and we define the order parameter to be (O2) = V2^ 2 \ following 

1(3 



the conventions of Ref. 9|. Although ip^ and ip^ do not have the same holographic 
interpretation in the deconstructed theory away from the continuum limit, our definitions 
provide an accurate measure of whether a nonvanishing profile of the ipj has been dynamically 
generated. 

IV. RESULTS 

In this section, we numerically solve the discretized theory of Sec. II III for N G 
{5, 10, 100, 1000}. In the case of iV = 1000, for the purpose of connecting with the con- 
tinuum theory, we fix the UV cutoff z\ — e = 10 _1 , while for other N we set e = a, the bulk 
lattice spacing. We fix the lattice spacing at the horizon an = 10 -5 . In all cases, we impose 
the following boundary conditions: 

0'i = -p = -l, ^ (1) = 0, <p N = 0, and V4r_i = ^-^iv , (4.1) 

where we have used the scaling symmetry of the theory to fix T c by setting p — 1 . Searching 
for nonvanishing ijjj solutions as we vary the temperature, we find there is a critical temper- 
ature for each N at which the operator (O2) condenses, as shown in Fig. [2j We note that 



there are numerous non-zero solutions for ipj but, as in the continuum model {q], we retain 
only the monotonic solutions. For N = 5, N = 10, N = 100, and N = 1000 we find the 
critical temperatures T c = 0.076 p 1/2 , T c = 0.088 p 1/2 , T c = 0.111 p 1/2 , and T c = 0.118 p 1/2 
respectively. Here we reintroduced the factor of p 1 / 2 as indicated by the scaling relations. It 
is interesting to see that for all N we retain the square root behavior of the phase transition. 
We find that the curves in Fig. [2]near T c are well fit by the form (0 2 ) = CnT^(1 — T /T c ) 1 / 2 . 
The values we obtain for the coefficients are Cinno = 127, C wo = 20, C w = 58, and C 5 = 93, 
as compared to the continuum value C = 144 [9]. 

We would also like to study the behavior of the complex conductivity away from the 
continuum limit. Using the expression for the current Eq. f)3.29p . we find the following 
expression for the conductivity: 

a _ J x _ f _ if i(Ax2 -A VVx )/a ^ 
E x A x wA UVa: 

As discussed in Sec. 1111} we impose the ingoing- wave boundary condition, Eq. ( I3.25p . n sites 
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FIG. 2: This figure displays the condensation phase transition at critical temperature T = T c . The 
top curve is the N = 1000 curve and below that the dashed curves from top to bottom around T c 
are N = 5, N = 10, and N = 100, respectively. The 1000-site curve includes a fixed UV cutoff e, 
which allows for a smooth continuum limit. 

away from the horizon. In particular, we set 

A xN - n = 1 and A xN _ n ^ = 1 - tUa , (4.3) 

JN-n-l 

where we have fixed the arbitrary normalization A x ^_ n = 1. Our choices for the shift n 
for N = 1000, 100, 10 and 5 are n = 20, 10, 2 and 2, respectively. With these boundary 
conditions for A X j we compute the real and imaginary parts of a, as shown in Figs. |3] and 
HI It is important to note that the graphs of Im[cr(a;)] in the superconducting phase exhibit 
a simple pole at u = 0. This implies the presence of a delta function contribution to the 
real part, Re[er(u;)] oc 5(oj), as discussed in Sec. [TTJ For iV = 10 there are a series of peaks in 
Re[cr(cj)] with a corresponding pole-like structure in the imaginary part as follows from the 
Kramers-Kronig relation Eq. (12. 17ft . Finally, the oscillatory behavior of Re[<r(c<;)] at large 
uj displayed in the N = 100 and N = 1000 plots is consistent with our expectation from 
Eq. (I2.19p given that our boundary condition is only approximately ingoing wave. 
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FIG. 3: Complex conductivities. The (N,T/T C ) values for the top and bottom rows are given 
by (1000,0.55) and (100,0.53), respectively. The dotted (red) curve is the normal (n) phase while 
the dashed (blue) curve is the superconducting (sc) phase. The vertical arrow represents a delta 
function. The solid (black) curve is the ratio: Re(<7 sc )/Re(<7 n ). 



V. CONCLUSIONS 



Dimensional deconstruction is a powerful technique for finding lower- dimensional theories 
that share some of the interesting phenomenological features of their higher-dimensional 
progenitors. In this paper, we have considered the dimensional deconstruction of four- 
dimensional holographic theories of superconductivity. We have shown how the AdS/CFT 
prescriptions for computing physical quantities of interest (for example, charge densities and 
currents) emerge in the lower-dimensional, discretized theory. We have also demonstrated 
that deconstructed theories with a relatively small number of lattice sites (here taken > 5) 
do indeed retain enough of the physics of the continuum limit to provide possible models of 
superconductivity on their own. 
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FIG. 4: Complex conductivities. The (N,T/T C ) values for the top and bottom rows are given 
by (10,0.54) and (5,0.50), respectively. The dotted (red) curve is the normal (n) phase while 
the dashed (blue) curve is the superconducting (sc) phase. The vertical arrow represents a delta 
function. The solid (black) curve (omitted for N = 10, for clarity) is the ratio: Re(<7 sc )/Re(cr n ). 

We view the results of this work as preliminary, in the sense that we have only considered 
the deconstruction of one of the simplest, s-wave holographic superconductors that were 
proposed in the first papers on the application of the AdS/CFT correspondence to this 
problem. Much of what we have found in the present work leads to questions and suggestions 
for future study: 

• Dimensional deconstruction also motivates new models of p- and c?-wave superconduc- 
tors. It would be interesting to study the phenomenology of those models. 

• In the present approach, the photon field is not included as a dynamical degree of 
freedom in the deconstructed theory. However, the theory can be easily modified by 
gauging the electromagnetic U(l) symmetry and adding an associated gauge kinetic 
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term. 

• The deconstructed model suggests the possibility of hidden local symmetries in the 
interactions between Cooper pairs and excitons. Perhaps absorption studies would be 
able to test the existence of these states and interactions. 

• If one deviates from a strict matching to the continuum holographic theory, then one 
has greater freedom in constructing a model that is based on a very small number of 
replicated gauge groups (for example, two or three). Would a two- or three-site model 
of superconductivity be phenomenologically viable? 

• The deconstructed model introduces temperature-dependent couplings inferred by a 
holographic model. It would be important to deduce these couplings instead from a 
microscopic description. 

• In the present work, the fluctuations of the link fields were ignored, consistent with 
the tacit assumption that these fields are heavy. However, these fields, which have no 
counterpart in the continuum holographic theory, might also have a physical interpre- 
tation. 

We look forward to considering these issues in more detail in future work. 
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Appendix A: A purely dynamical ip 

To simplify our discussion, it was assumed earlier that the ipj fields at the ends of the 
moose, ipi and ipw, were background fields. Here we show that the same boundary conditions 
would be obtained in the continuum limit if we allow all the ipj to be dynamical fields. The 
Lagrangian that is relevant in this case is a minor modification of Eq. ( 13. 14f> : the third sum 
in this expression, which runs from j = 2 to N — 1, is changed to one which runs from j = 1 
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to N (with <2tv = clh) for the first two terms in square brackets and from j = 2 to N for the 
mass squared term. We also shift Zn away from its previous value by e#, which serves as a 
cutoff to regulate the divergence in 1/fpj. Eq. ( I3.16P remains unaltered, but now there are 
two additional equations of motion, for ipi and 



a Ji 







and 



—^N-l 

a H 



"N-l 



N 



N 







(Al) 



(A2) 



where V'jv-i = (4>n — ipN~i)/ a H- These equations can be viewed as dictating the boundary 
conditions for Eq. (I3.14p . which, in our previous approach, were chosen freely to mimic the 
boundary conditions of the continuum theory. Eq. (IA1I) reduces in the continuum limit to 
the boundary condition ip[ = 0. Since 4>n = 0, the second term in Eq. flA2j) vanishes (1//jv 
is finite for finite cutoff e#). Hence, 

2anz%_ 1 



N-l 



N 



(A3) 



In-izJ? _ 

The quantity in brackets can be expanded in powers of an while maintaining C a^; one 
finds 



N-l 



3z H 



34 



a H + 0(a 2 H ) 



N 



(A4) 



In the limit an — > we recover the boundary condition in Eq. (13.191) . In summary, the effect 
of treating the ipj as dynamical fields everywhere is to modify the form of the boundary 
conditions away from those of the continuum limit: 

2 2 
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